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Abstract 

To describe the flow of a miscible quantity on a network, we introduce 
the graph wave equation where the standard continuous Laplacian is re- 
placed by the graph Laplacian. This is a natural description of an array 
of inductances and capacities, of fluid flow in a network of ducts and of 
a system of masses and springs. The structure of the graph influences 
strongly the dynamics which is naturally described using the basis of the 
eigenvectors. In particular, we show that if two outer nodes are connected 
to a common third node with the same coupling, then this coupling is an 
eigenvalue of the Laplacian. Assuming the graph is forced and damped 
at specific nodes, we derive the amplitude equations. These are analyzed 
for two simple non trivial networks: a tree and a graph with a cycle. 
Forcing the network at a resonant frequency reveals that damping can be 
ineffective if applied to the wrong node, leading to a disastrous resonance 
and destruction of the network. These results could be useful for complex 
physical networks and engineering networks like power grids. 



1 Introduction 



The flow of a scalar quantity in a network is an important problem for funda- 
mental science and engineering applications. The latter include gas, water or 
power distribution networks. Other examples are a simplified version of road 
traffic and the flow of nerve impulse in the brain. The static aspect of the prob- 
lem has long been studied within the framework of operations research, see for 
example [T]. However in many cases the dynamic character is crucial. Take for 
example the prediction of traffic jams in a given road network or the prediction 
of a flood in a river basin. 

The basic equations describing the flow of a miscible quantity are the well- 
known conservation laws of mathematical physics. These laws are building 
blocks for studying physical systems. They are universal and can be found in 
mechanics, in electromagnetism .. Each conservation law has the form of a flux 
relation 

ut + d^q = 0, (1) 

where the flrst term is the time derivative of a density u and the second one is the 
space derivative of the flux q along the relevant coordinate. The most important 
conserved quantities in mechanics are the mass, momentum and energy. In 
electromagnetism the current and voltage obey conservation laws. 

To describe the flow of a given quantity on a network it is natural to try and 
generalize these conservation laws. For that we introduce the graph represent- 
ing the network and the generalized gradient V or its transpose, the incidence 
matrix. To write this it is important to orient the branches of the graph in 
a flxed way. This can be arbitrary. An important class of flows are the ones 
such that there is no dissipation along the branches but only at the vertices. A 
typical example is a small power grid for which the power line dissipation can 
be neglected and where the only power input and outputs occur a given nodes. 
For some models it is possible to reduce the dynamics to what has been called 
a graph wave equation by Friedman and Tillich [5] . Here the usual Laplacian is 
replaced by its discrete analog the graph laplacian V-'^V. 

Here we show how the conservation laws lead to the graph wave equation. We 
illustrate this in three different physical contexts: a network of inductances and 
capacities, its equivalent mechanical analog represented by masses and springs 
and an array of fluid ducts. For the latter, note the study by Maas 3 who 
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considered graphs obtained by linking elementary graphs. He established in 
particular inequalities for the first non zero eigenvalue of the Laplacian of these 
graphs. Note also that the static problem of gas or fluid transport is usually 
solved using optimisation techniques [4]. The graph laplacian was also consid- 
ered recently by Burioni et al for the thermodynamics of the Hubbard model to 
describe static configurations of a Josephson junction array [S]. Another prob- 
lem considered by that group is to use the graph as a controlled obstacle to 
obtain a desired reflection of discrete nonlinear Schrodinger solitons [6]. Here 
we adopt a different point of view. We consider that the network is fixed and 
is submitted to forcing and damping on specific nodes. This formulation is par- 
ticularly interesting because the graph Laplacian i.e. the spatial part of the 
equation is a symmetric matrix so that its eigenvalues are real and its eigenvec- 
tors are orthogonal. It is then natural to describe the dynamics of the network 
by projecting it on the basis of the eigenvectors. Here we examine these normal 
modes for two simple graphs, a tree and a graph with a cycle. We write down 
amplitude equations for the normal modes of the system. Finally we force the 
graph on a given node and damp it on another. We choose the forcing frequency 
so that the system resonates. We illustrate two specific effects: first that damp- 
ing can be ineffective if applied to the wrong node. The second is that we can 
have a multiple eigenvalue corresponding to different eigenvectors. Then de- 
pending how the system is excited we can get a different response. We consider 
a periodic forcing for simplicity. The result for other types of sollicitations can 
be derived from this study using a superposition argument because the system 
is linear. 

The article is organized as follows, section 2 presents the derivation of the graph 
wave equation in different physical contexts. In section 3 we compute the eigen- 
modes for two specific simple graphs, a tree and a graph with a cycle. Section 
4 introduces the equations for the amplitudes of the normal modes when the 
graph is forced. These equations are analyzed in section 5. There we force the 
simple graphs of section 3 at resonance and analyze their response. We conclude 
in section 6. 
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Figure 1: Schematic drawing of the complex graph oscillator. The arrows on 
the branches are oriented arbitrarily. 

2 The model: graph wave equation 

We now introduce the basic notions from graph theory following the presentation 
of 0. A graph Q{V,E) is the association of a vertex set V and an edge set E 
where an edge is an unordered pair of distinct vertices. We assume the vertices 
and edges to be numbered V — {1, 2, . . .n} and E = {1, 2, . . . m}. The latter 
are oriented with an arbitrary but fixed orientation. We consider for simplicity 
only simple graphs which do not have multiple edges. Fig. [Ijshows such a graph 
with n = A vertices and m = 3 edges. The basic tool for expressing a flux is the 
so-called incidence matrix C(n,m) defined as 

Cxe = — 1 edge e starts from vertex x, (2) 
Cxe = 1 edge e finishes at vertex x, (3) 
Cxe = otherwise (4) 

For the example shown in Fig. [TJ we have 
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The transpose C'^ = V is a discrete differential operator ( gradient of graph). 
To see this consider a function f : V R. The vector (V/)(e) is the difference 
of the values of / at the end points of vertex e (with orientation). In the example 
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above, we have 



-1 1 





V 



//A 

h 



(h 




h 
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which is the discrete gradient of / associated to the graph. 

We now consider the specific inductance-capacity electrical network shown 
in the left panel of Fig. [2l The equations of motion in terms of the (node) 
voltages and (branch) currents are the conservation of current and voltage 



Cvt - 
Lit 



- Vw = 0, 
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are respectively the diagonal matrix of capacities and the diagonal matrix of 
inductances and s represents the currents applied to each vertex, (similar to 
boundary conditions in the continuum case). Note that the equations © also 
describe in the linear limit, shallow water waves in a network of canals and fluid 
flow in a pipe network. In the first situation v corresponds to a surface elevation 
and g to a flow. For the second example v is the pressure and q the flow. From 
the two equations ^ one obtains the generalized wave equation. For this take 
the derivative of the 1st equation and substitute the second to obtain 
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A similar equation arises when describing the other physical system shown 
in the right panel of Fig. [21 the collection of four masses , i = 1 — 4 connected 
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Figure 2: Two different physical networks corresponding to the graph shown in 
Fig. [T] The left panel shows a network of inductances and capacities and the 
right panel shows the mechanical analog in terms of masses and springs 

by springs of stiffnesses ki, i = 1 — 3. Here the variables are the displacements 
Hi, i = 1 — 4 of each mass. The equations of motion are 

miy'i = ki {y2 - yi) 

TO22/2 = ki (yi - y2) + k2 (j/3 - 2/2) + ^3 (2/4 - 2/2) 
msy's ^ k2 {y2 - ys) 
^ TO42/4 ^ ks {y2 - yi) ■ 

Notice the correspondence capacities / masses and inverse inductances / stiff- 
nesses. The matrix on the right hand side is symmetric. This symmetry has 
important consequences as shown below. 

In the following we will consider the graph shown in Fig. [3] which includes 
an additional branch between nodes 3 and 4. We assume that the masses are 
equal, that fci = fca = fc, k2 = a^k and k^ = j5'^k. 
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Figure 3: The specific network that will be studied throughout the article. 



We normalize times by the natural frequency 



— , tl = ut. 

m 



(10) 



Omitting the primes, the equations can be written in matrix form as 
( yi\ ( -I 1 \ / yi \ 



2/2 _ 1 - 
2/3 ~ 
V 2/4 / V 

which we will write formally as 



2 - a a 1 
a —a — j3 l3 
1 (3 -1-/3/ 



Ytt = GY, 



2/2 
2/3 
V 2/4 / 



(11) 



(12) 



where G is the graph Laplacian [7], the equivalent of ([9]). The matrix G is 
symmetric. If we had assumed different masses on the nodes, we would have 
lost this property. For electrical networks this corresponds to the same capacity. 
In the rest of the article we keep the masses the same. 

Notice that G is a singular matrix since the sum of its lines (resp. columns) 
gives a line (resp. column). Therefore G will have a zero eigenvalue which 
corresponds to the Goldstone mode. The solutions of the linear system GY = S 
are given up to a constant. The linear evolution problem (jl2l) gives rise to 
periodic solutions Y(t) = Zexpiujt where the Z verify the spectral problem 



GZ 



2 ry- 

-U! Z. 



(13) 



Since the matrix G is symmetric, its eigenvalues are real and the eigenvectors 
are orthogonal. They provide a basis of i?" which is adapted to describe the 



7 



evolution of Y on the graph. Specificahy we arrange the eigenvalues = — 
as 

|Ai| =0< IA2I < ••■ < |A„|. (14) 
We label the associated eigenvectors vi,V2, ■ ■ ■ ,Vn- These verify 

Gvt = -Lofvi. (15) 

and are orthogonal with respect to the standard scalar product. We then nor- 
malize them so < Vi^Vj >= Si^ where Sij is the Kronecker symbol. 

It is natural to write the equation of motion (jl2p in terms of the amplitudes 
of the normalized eigenvectors Vi, Y — J2i ctiVi- We obtain the standard result 

di + u'fat = 0, (16) 

so that the normal modes do not exchange energy. The problem is then Hamil- 
tonian with 

iJ=i^d^+a.,W. (17) 

i 

When a perturbation acts on the system the equations need to be modified. 
This is the object of section 4. 

In this article, to make things precise, we have chosen to study the network 
shown in Fig. [3] The dynamics on this simple yet not trivial example already 
shows some specific effects like the influence of damping and the multiplicity of 
eigenvectors for a given eigenvalue. When /3 — the graph is called a tree (a 
connected graph with no cycle). We will consider this situation first. We will 
see that for a graph with a cycle (/3 ^ 0) the dynamics will be quite different 
than for the tree. 

3 Computation of the eigenmodes 
3.1 The case of a tree 

Throughout this section we assume that the branch 4 is absent so that /3 = 0. 
Let us first analyze the degenerate case a — 1. Then the eigenvalues are 

A4 = -4 A3 = -1 A2 = -1 Ai = 0. (18) 
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As expected we have the Goldstone mode Ai = with its fiat eigenvector 
(1, 1, 1, 1). The eigenvector associated to A3 = A2 are (1, 0, 0, —1) and (0, 0, 1, —1). 
The former is the antisymmetric mode which will be preserved when perturbing 
the graph around y^. The eigenvector for A4 = —4 is (1, —3, 1, 1). The eigen- 
value -1 is double. The case a = 1 is a special situation unlikely to occur for 
real systems for which one branch will always have a different stiffness from the 
other. For this reason we will assume a ^ 1. 

For a general a the eigenvalues and eigenvectors can be computed analyti- 
cally. As previously we find the Goldstone mode Ai = — cj^ = with the usual 
eigenvector. The other eigenvalues and eigenvectors are 

A4 = = -i (V4a2 _4q, + 9 + 2a + 3) , (19) 
vl = {^^~\ (V4a2 -4a + 9 + 2a + l) , ^ (V4a2 -4a + 9 + 2a - 3) , 1^ , (20) 
and 

A3 = ^ i (^V4a2 -4a + 9 - 2a - 3^ , (21) 
vl = ^1, i (\/4a2 -4a-h9 - 2a - l) , (V4a2 -4a + 9 - 2a + 3) , 1^ , (22) 
and 

A2 = -ujI = -1, (23) 
«2^ = (1,0,0,-1). (24) 

The frequencies uJi are plotted as a function of the parameter a in Fig. Note 
how the eigenvalue W2 is independant of a, uj2 = 1- We will explain this lower 
in the section. 
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Figure 4: Plot of the eigenvalues w,;, j = 1 — 4 as a function of a for the graph 
shown in Fig. 3 with (3 = i.e. a tree. 

Notice how the mode V2 is unchanged, it is shown schematically in Fig. [SJ 



Figure 5: Schematic representation of the constant mode U2, corresponding to 

Fig. ini shows the plot of the components of the non trivial eigenvectors as 
a function of a. Note how the nodes 1,2 and 4 are in phase for ^3 for a < 1. 
For a > 1 node 2 becomes out of phase with nodes 1 and 4. Interestingly for 
a = 1 the second component of W3 = 0, this will have important consequences 
for the dynamics as we will see below. The eigenvector V4 (right panel of Fig. 
[6]) is simpler, nodes 1,3 and 4 are always in phase. The modes V2 and ?;3 evolve 
continuously from a = 0.1 to 0.9. Only the amplitudes change. When a = 5 
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the mode V2 has changed, the nodes 2 and 3 are in phase while the nodes 1 and 
4 are out of phase. 




Figure 6: Plot of the components of the non trivial eigenvectors (left panel) 
and W4 (right panel) as a function of a. 
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Figure 7: Schematic representation of the non trivial eigenvectors W3,W4. The 
top panel is for a = 0.1. It shows the modes V3 corresponding to W3 = 0.36 (left) 
and V4 corresponding to UI4 = 1.75 (right). The middle panel is for a = 0.9 and 
shows the modes vs corresponding to ^3 = 0.964 (left) and V4 corresponding 
to LJ4 = 1.9671 (right). The bottom panel corresponds to a = 5. It shows the 
modes vs corresponding to W3 = 1.335^ (left) and V4 (right) corresponding to 
u)4 = 3.35. 



The Figs. [S] and [7] show a schematic representation of the non trivial eigen- 
vectors for a = 0.1, 0.9 and 5. 

3.2 A general result for the eigenvalue 1 

The tree studied in the previous section has a Laplacian whose eigenvalue is 
equal to 1 corresponding to an eigenvector independent of a. This is a general 
property for any graph such that two outer nodes (termed leaves in graph theory) 
are connected to a common node with the same coupling. Let us assume this 
coupling is 1 for now. 

We consider the situation shown in Fig. [8] where the outer nodes 1 and 2 
are connected to a common node 3 which has p connections to the rest Q' of the 
graph. 




Figure 8: A graph such that two outer nodes are connected to a common node 
with the same coupling. 



The characteristic polynomial associated to the Laplacian matrix G is 



det(G - XI) 



1-A -1 

1-A -1 

-1 -1 P + 2-X -1 ... -1 
0-1 



-1 








(25) 
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Adding all the lines to the first line, we obtain —A for all coefficients of the line. 
Factoring it, we get 



det(G - XI) = -X 



1 1 

1 - A 
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-1 



-1 




-1 P + 2-X -1 
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We then add the first line to the third line to obtain 
11 11 
1-A -1 

J3 + 2-A + 1 
-1 

det(G - XI) ^ -A 



1 





We can then expand using the first column and the second column successively 
so that the factor (1 — A) appears in det(G — XI). 

Remark that if the couplings between nodes 1 and 3 and nodes 1 and 2 are 
equal to a instead of 1, we would get a as an eigenvalue. 

Using a similar argument it can be shown that if k leaves are connected to a 
common node with the same coupling a then a is an eigenvalue of multiplicity 
yfc-l. 

The eigenvector associated to the eigenvalue 1 (or a) is such that 3:3 = 0. 
This can be seen by inspection of the first line of (G — XI), see formula (P5|) . 
The third line of (G — XI) shows that 



ier(3) 



Xi = 0, 



where r(3) is the set of the nodes adjacent to node 3. 
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3.3 The graph with a cycle 

We now add branch 4 to the graph and form a cycle. We assume for simphcity 
a = /3. For this more complex graph, the eigenvalues and eigenvectors must 
be computed numerically. We have used Matlab. As for the tree case we have 
chosen a = 0.1, 0.9 and 5. The plot of the eigenvalues as a function of a is 
shown in Fig. [HI 

4 
3 

S 2 
1 


1 2 3 4 5 

a 

Figure 9: Plot of the eigenvalues w,;, j = 1 — 4 as a function of a for the graph 
shown in Fig. 3 with a = /3 i.e. a graph with a cycle. 

As expected from the theory 8 , the eigenvalues of the graph with a cycle ^ 
and the eigenvalues of the tree A are interlaced such that 

= Ai = ^1 < A2 < /i2 < A3 < < A4 < TOU4. (26) 

Note that A4 = 114 for a < 1. The fact that A2 and A3 are close confines ^2- 
This has important consequences for the dynamics. 

The dependency of the eigenvectors on the coupling parameter a is shown 
in Fig. [TOl In the next section, we study it in detail and compare it to the one 
for the tree. 
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Figure 10: Plot of the components of the non trivial eigenvectors W2, wa, i'4 from 
top to bottom as a function of a. 



3.4 Comparison of the eigenmodes between the tree and 
the graph 

Let us compare the modes for different vaiues of tlie coupling parameter a. 
Consider tire small value a = 0.1. For the mode 2, The tree and the graph 
with a cycle show similar eigenmodes 2 and 3. For the 4th eigenmode the cycle 
breaks the symmetry between the nodes 1 and 4 but this change remains small 
as shown in table 4 . 




Figure 11: Strong coupling a = 0.9. Comparison of the eigenmodes 2 , 3 and 4 
from left to right for the tree (top panels) and the graph with a cycle (bottom 
panels). 

The modes 2 3 and 4 for a = 0.9 are shown in Fig. [11] for the tree (top 
panels) and for the graph with a cycle (bottom panels). For this value of a the 
branches 2 and 4 of the graph with a cycle are almost equivalent so that we 
expect a large difference between the modes for the graph with a cycle and the 
modes for the tree. This is indeed the case. The modes 3 are close as shown 
by the frequencies. On the contrary the mode 2 for the graph with a cycle only 
involves nodes 3 and 4, its frequency is almost double the one for the mode 2 for 
the tree. The mode 4 for the graph with a cycle couples nodes 3 and 4 together 
in opposition to node 1 and its frequency is almost equal to the one for the tree. 

Finally, we discuss the high value of the coupling a = 5. Fig. [TH shows 
the modes 2, 3 and 4 for the tree (top panels) and for the graph with a cycle 
(bottom panels). The mode 2 has a frequency of 1.3 for the tree and 3.90 for 
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the graph with cycle. There the hnk between nodes 3 and 4 is inactive. For the 
mode 3 activity switches from the node 3 to the node 4 as one adds the Hnk 
between 3 and 4. Finally note how the mode 4 is modified for the graph with a 
cycle. The nodes 3 and 4 are now strongly coupled and oscillate in phase. The 
frequency 1.10 remains close to 1. 




Figure 12: a = 5, comparison of the eigenmodes 2 , 3 and 4 from left to right 
for the tree (top panels) and the graph with a cycle (bottom panels) . 



4 Forcing the network : amplitude equations 

For a general (nonlinear) evolution problem of the form 

Ytt^GY + N{Y), (27) 

it is natural to expand Y using the eigenvectors as 

Y(t) = ai(t)vi + a2{t)v2 + ■■■ + an{t)vn. (28) 

Inserting (pS)) into ([?7|) and projecting on each mode Vi we get the system of 
coupled equations 

aitt+ujlai = <N{Y)vi>, (29) 
a2tt+^la2 = <N{Y)v2>, (30) 

(31) 

Untt + ^lan = < N{Y)Vn > . (32) 
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We expect this decomposition to be more adapted to describe the dynamics of Y 
on the graph. In particular it should explain some of the unexpected couplings 
that are observed between the modes. 

Let us now assume that the network is forced at some node n f and damped 
at some node Ud- The motion can be represented as 

Ytt = GY + Md{~dYt) + Mffl, (33) 

where the (71,71) matrices Md, Mf are everywhere except for Md{nd,nd) = 
1, Mf{nf,nf) = 1. The vector 1 is 1 = (1,1,...,!)^. Assuming the linear 
combination for Y ([25)1 and projecting we get 

n 

dj = -uj^aj -d^< MdVk\vj > cik + f < Mfl\vj > . (34) 
fe=i 

In terms of components we obtain 

n 

a J = -w|aj - dv^'' ^ Vk''°-k + /^^"^ ; J = 1, " (35) 
fc=i 

where vj'' , vj^ are respectively the 71^, Uf components of the normal vector vj. 
From these equations one can see that exciting one node will cause disturbances 
to propagate all through the network in a precise way. The forcing will act on 
mode j only if the component vj^ ^ 0. The damping of mode j will be effective 
only if v""^ ^ 0. In the next section we will show examples of the dynamics 
of the network when it is excited periodically at a given node and damped at 
another. We will examine this for both a tree /? = and a graph with a closed 
cycle. 

5 Numerical results: forcing the network 

As an example of the dynamics of the network, we consider that the graph is 
forced periodically at a given node ?i/ and damped at a node ti^, this happening 
for a time duration 100 < t < 300. This could model an electrical power 
grid where the power stations input energy. The nodes where damping occurs 
correspond to cities where the energy is absorbed. 

As seen above the coordinates of the eigenvectors will determine whether a 
given node will contribute to the mode amplitude or not. So a first interesting 
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effect is that damping can be completely ineffective if it is applied to the wrong 
node. To illustrate this effect we choose a tree with a = 0.1, (3 — and the 
same graph but with a cycle so that a = 0.1, (3 — 0.1. The first set of figures 
correspond to exciting the tree on node 4 and damping it on node 1 and on 
node 2. For the damping on node 1 the corresponding amplitude equations are 

(36) 
(37) 
(38) 
(39) 

The 2nd equation is resonant if / = sint {uj — I ) but it is damped. This 
gives rise to an increase of the amplitude of mode 2 with a saturation as shown 
in the left panel of Fig. [T31 The same happens if the node 4 is damped. The 
final value of 02 can be obtained by a Green's function approach on the 4th 
differential equation. This equation can be written 

o + w^o = f{t) - da. (40) 

The Green's function G{t — to) solves 

G + ujlG + dG^S{t-to), 

where S{t — t^) is the Dirac function at to. From ^ we have 

sin x/wn - ^(t - to) 

A ,2 

y^o - — 

The solution is the given by 

a(t)= / G{t-tQ)f{to)dto. (42) 
^0 



On the contrary when damping is applied to nodes 3 or 4 the equation for 02 
becomes 

^2 + 02 = -0.7/, (43) 

which is not damped at all. The amplitude of the mode 2 then grows linearly 
as expected. This is shown in the right panel of Fig. [T31 The final amplitude in 
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Figure 13: Example of the influence of damping applied to different nodes, the 
graph is a tree with a = 0.1, /3 = 0. Plot of 02 (t) (blue online), a-^ii:) (red 
online) and 04 (i) (green online) when the node 4 is excited with a frequency 
10 — \. The left (resp. right) panel corresponds to the node 1 (resp. 2) being 
damped. 

Fig. [13] can be easily found by seeing that the solution of the linear resonance 
equation 

a + a = / sint, (44) 

is 

a{t) = -^-tco%t. (45) 

In the example above one gets 

maa;(a2(300)) = -^(300 - 100) = 70, 

which is in excellent agreement with the value on the right panel of Fig. 1131 

When the graph has a cycle so that /3 = 0.1 we get a slightly different 
picture because the linear resonance observed previously is not exact. There 
is a small damping due to the non zero second coordinate of the eigenvector 
v-j,. We excite node 4 at a frequency uj = \ and damp the nodes 1,2,3 and 4 
respectively. We find two groups of behaviors shown in Fig. [TH the left panel 
corresponds to damping on nodes 2 and 3. One can see the typical beat at 
frequency (^3 — a;)/2 which yields a half-period T/2 = 157. The mode 2 which 
is the initial condition is strongly reduced. When the damping is applied to the 
nodes 1 or 4 we have a damped resonance. This is shown in the right panel of 
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Fig. [TH Note that the amplitude of mode 3 is practically constant during and 
after the forcing/damping region. 
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Figure 14: Similar plot as in Fig. [T3] except that the graph now has a cycle 
a — 0.1, P = 0.1. Plot of a3{t) (blue online) and a2{t) (red online) when 
the node 4 is excited with a frequency cj = 1. The left (resp. right) panel 
corresponds to the node 2 (resp. 1) being damped. 

Finally note that if we excite nodes 2 or 3 with uj = 1 we do not get any 
significant response. This is because V2 has zero components on nodes 2 and 3. 

Another effect occurs because close frequencies correspond to two different 
eigenmodes. So depending on the way the network is excited or damped we get 
the two eigenmodes or just one. As an example consider the tree graph with 
a = 0.9, (3 — 0. As shown in Fig. |4]we have two close eigenvalues 
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Figure 15: Plot of 02 (i) (blue online) and 03 (t) (red online) when the node 4 is 
excited with a frequency uj = 1. The panels correspond respectively to damping 
the 1st node (left), damping the second node (middle) and damping the 3rd 
node (right). 

When node 1 is damped, both modes V2 and U3 can be excited because they 
have non zero components on that node. This is shown in the left panel of 
Fig. [m On the contrary when node 2 is damped (middle panel of Fig. [T5|) no 
damping occurs for the amplitude 02 causing an unbounded linear growth. The 
mode V3 is excited but in a much smaller way because it is weakly damped, since 
W3 has a small non zero component on node 2. When node 3 is damped, the 
mode 3 is strongly damped so that there is practically only mode 2. Therefore 
one sees that applying damping to node 3 will result in mode 2 only being 
excited while applying damping on node 1 will result in the presence of both 
modes 2 and 3. 

6 Conclusion 

To describe the flow of a miscible quantity on a network, we introduced the 
graph wave equation where the standard continuous Laplacian is replaced by 
the graph Laplacian. We showed that a natural example is an electrical network 
of inductances on the branches and capacities on the nodes. This system can 
also describe shallow water waves on a network of canals or fluid flow in a 
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network of pipes. There is also a mechanical analog in terms of masses and 
springs. 

Since the graph Laplacian is a symmetric matrix, its eigenvectors are orthog- 
onal and provide a natural basis to describe the flow in terms of amplitudes on 
each mode. We derived such amplitude equations when the network is forced 
and damped on a given node. The eigenvalues and components of the eigenvec- 
tors are important elements of these amplitude equations. We analyzed them 
for two simple non trivial networks: a tree and a graph with a cycle. For a tree 
such that at least two outer modes are connected to a common third node with 
the same coupling a , we get the general result that one eigenvalue is equal to 
a. 

The numerical analysis of the amplitude equations shows in particular that 
damping can be ineffective if applied to the wrong node. This could cause 
disastrous resonance and destruction of the network. Another effect is that 
we have multiple eigenvalues corresponding to different eigenvectors. Therefore 
exciting the system on different nodes can cause one or the other eigenvector to 
appear. 

These results could be useful for complex physical networks like arrays of 
Josephson junctions. They could also have important applications in engineering 
situations like for a power grid. 
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